Metallic phase of the quantum Hall effect in four-dimensional space 
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We study the phase diagram of the quantum Hall effect in four-dimensional (4D) space. Unlike 
in 2D, in 4D there exists a metallic as well as an insulating phase, depending on the disorder 
strength. The critical exponent v « 1.2 of the diverging localization length at the quantum Hall 
insulator-to-metal transition differs from the semiclassical value v = 1 of 4D Anderson transitions 
in the presence of time-reversal symmetry. Our numerical analysis is based on a mapping of the 4D 
Hamiltonian onto a ID dynamical system, providing a route towards the experimental realization 
of the 4D quantum Hall effect. 
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To understand our three-dimensional world it is help- 
ful and rewarding to ask how physics would differ in 
other dimensions. The interest in two-dimensional and 
four-dimensional worlds goes back to the 19th century 
[IHH , inspired by the development of non- Euclidean ge- 
ometry. In modern times, nanotechnology has brought 
us the 2D "Flatland" in semiconductor heterostructures 
and in graphene [H, while ultracold atomic lattices offer 
a way to simulate 4D dynamics by encoding the extra 
dimension in an internal degree of freedom of the atoms 

These studies of physics in other dimensions are par- 
ticularly interesting if they reveal effects that lack a 3D 
counterpart. The integer quantum Hall effect is a cele- 
brated example, first discovered in 2D Q, then general- 
ized to 4D and 8D 7, 8:], and now known to exist in any 
eiien-dimensional space 0-[H|- In each case an integer- 
valued topological invariant (a Chern number C [1Ch12| |) 
counts the number of gapless surface modes, observable 
as the quantized Hall conductance in 2D. Disorder lo- 
calizes single-particle excitations in the bulk but not on 
the surface, producing an insulating quantum Hall fluid 
bounded by a conducting surface. 

Topological phase transitions (at which C switches 
from one integer to another) have been extensively stud- 
ied in the 2D quantum Hall effect [l3j], but not for the 
higher dimensional generalizations. As we will show here, 
in 4D the quantum Hall insulators with a different topo- 
logical invariant are separated in phase space by a metal- 
lic region, so that a topological phase transition is com- 
posed of a pair of metal-insulator transitions. In 2D an 
intermediate metallic phase requires the presence of time- 
reversal symmetry (as in the quantum spin-Hall effect 
[li]]), but in 4D there is no such requirement. We find 
that the 4D quantum Hall insulator-to-metal transition 
has a different critical behavior than the 4D Anderson 
transitions in time-reversally-symmetric systems [TH, . 

To determine the phase diagram of the 4D quantum 
Hall effect we proceed as follows: We start from a Dirac 
Hamiltonian 0, [To| of electrons with 4D momentum p 
coupled to SU(2) spin and valley degrees of freedom cr, 



t. This model has time-reversal symmetry, so it is in the 
universality class of the quantum spin-Hall effect (class 
All in the classification of Anderson transitions [ll]). In 
order to bring it in the universality class A of the quan- 
tum Hall effect, we add time- reversal-symmetry breaking 
disorder. Following the lines set out in Ref. for the 2D 
quantum Hall effect, we then map the 4D Hamiltonian 
problem onto a ID dynamical system. Here we use this 
dimensional reduction to arrive at an efficient numerical 
simulation, but in principle this method could be used to 
simulate the 4D quantum Hall effect in ultracold atomic 
lattices [2(J. 

The Dirac Hamiltonian has the form [10, \2M 



2 arctan \v\ ^ 
H(p) = — > ViT. 
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(1) 



v = K ^fj, + ^2 c ospi, sinpi, sinp 2 , sinp 3 , sinp 4 J , (2) 
T = (t x ® cr , r 2 <g) er , i~ y ® cr x , r y ® a y , r y ® a z ). (3) 

The vectors v, T have components Vi,Ti (i = Q, 1, 2, 3, 4). 
The Pauli matrices a ai r a (a = x, y, z), with unit matri- 
ces <To, To, act on the spin and valley degrees of freedom. 
The two parameters K, \i are real numbers. 

The Hamiltonian (JXJ) is a 4 x 4 matrix dependent on 
the momentum vector p — (jp\ ,P2,P3,Pi), in the Brillouin 
zone \pi\ < 7r of a 4D square lattice (lattice constant 
and h set equal to unity). It satisfies the time- reversal 
symmetry relation 

H(p)=TH(-p)T-\ T = (i<T y 8 t x )K, (4) 

with JC the operator of complex conjugation. The spec- 
trum consists of a pair of twofold degenerate bands 



E±(p) = ±2 arctan \v\, 



(5) 



with a band gap centered at the Fermi level E = 0. The 
topological invariant (the second Chern number) depends 
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on the effective mass parameter //, according to [T(l [HJ 

[3 sign if 0< | M | <2, 
C= ^-sign(/x) if 2 < <4, (6) 
[0 if | M | >4. 

Changes in C are signaled by a closing and reopening 
of the band gap at time-reversally invariant momenta 
[pi = —pi (mod 27r)], where v(p) vanishes. 

We now include a spin-dependent potential V(x), with 
position operator x = (xi,X2,x 3 ,x 4 ) = id /dp, to break 
both translational invariance and time-reversal symme- 
try. To arrive at a dynamical system the potential is 
added stroboscopically in time, 

00 

H(t) = V(x)+H(p) ]T 5(t-n). (7) 

n— — 00 

The time evolution of a wave function over one period is 
given by the Floquet operator 

jr = e -iH( P ) e -iV(x)_ ( g ) 

(We measure position in units of the lattice constant and 
set the stroboscopic period equal to unity.) Periodically 
time-dependent models of this type are called quantum 
kicked rotators [23| . 

A first advantage of the stroboscopic description is that 
a simple separable potential, 

4 

V(x) = Vi(xi) - ^uiiXi, V 1 (x 1 )=x 2 1 T z ®a , (9) 

i=2 

can be used to simulate d-dimensional disorder — pro- 
vided the frequencies uji/2ir are incommensurate irra- 
tional numbers [U, l25j . Note that the t z matrix anti- 
commutes with T from Eq. so this potential breaks 
time-reversal symmetry. 

A second advantage is that the <i-dimensional dynam- 
ics with a single stroboscopic period can be mapped ex- 
actly onto a 1-dimensional dynamics with d stroboscopic 
periods [26j . For that purpose we choose an initial wave 
function of the form 

4 

*(p, t = 0) = V>(P! , t = 0) JJ S(jh - a,). (10) 

During one period the momentum pi (i = 2,3,4) is in- 
cremented to pi + Ui (mod 27r), so we may replace the 
4D dynamics by a ID dynamics with a time-dependent 
Floquet operator 

J~(t) = e ~ iH (P 1 > UJ2t + a2 ' UJ: > t + a3 ' UJ4t + a < L ) e - iV i( x i) (11) 

We introduce a Bloch number q and calculate the time 
dependence of the wave function ijj(pi,t) — e iqpi x(pi,t). 
The state x(j>l> t) is a 27r-periodic function of pi, so it is a 
superposition of a discrete set of eigenstates e" npi of x\. 
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FIG. 1: Phase diagram of the 4D quantum Hall effect (broken 
time-reversal symmetry, class A). The colour scale represents 
T> — A 2 jsft at t = 10 5 . The phase boundary separating metal 
and insulator appears approximately at T> « 0.08 (white). A 
more precise determination of the phase boundary from the 
scaling analysis is indicated by the red curves (with circles 
marking the two transitions analyzed in Table d}. Because of 
the ±fj, symmetry only positive values of fj, are shown. 

We truncate the set to M states, m £ {1, 2, . . . M }, with 
periodic boundary conditions at the end points. The cal- 
culation of ip at integer t then amounts to subsequent 
multiplications of the Floquet operator, which can be 
done efficiently using the Fast Fourier Transform algo- 
rithm Ranges of M up to 10 3 and t up to 6 • 10 7 are 
included in the simulations. 

The metallic and insulating regions in the phase dia- 
gram are obtained by initializing x( x ii a ^ %i — 3- n d 
calculating the mean square displacement at time t, 



A 2 W =^(m-m ) 2 |x(^i =m,i)| 2 . (12) 

n 1 

The overbar indicates an average over initial conditions 
and over the Bloch number q. In the large-t limit one has 
A(t) — > constant in the insulating phase (localization) 
and A 2 (i) oc t in the metallic phase (regular diffusion). 
On the critical line separating these two phases A 2 (t) oc 
y/i (anomalous diffusion). The resulting phase diagram 
as a function of effective mass /x and disorder strength 
1 /K is shown in Fig. [TJ 

The clean limit of the model is reached for K — > 00 
[because then V(x) is negligibly small relative to H(p)]. 
In this limit the system switches between topologically 
distinct insulating phases at = 0,2,4, see Eq. ((6J. 
Disorder introduces a metallic phase in between the in- 
sulating phases, so that now a change in C happens via 
two metal-insulator transitions. At fixed /i the metallic 
phase exists for a finite range of disorder strengths, van- 
ishing both for weak and strong disorder. The insulat- 
ing phase that appears at the largest disorder strengths 
(1/K > 10) is disconnected from the insulating phases 
at smaller disorder, but presumably it is topologically 
trivial (C = 0). 
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FIG. 2: Finite-time scaling of the mean-square displacement 
at the quantum Hall insulator-to-metal transition (left circle 
in Fig. [T]). The data points result from the numerical simula- 
tion, the curves are fits to the scaling law (|13|) . The crossing 
identifies the critical point fi c , separating the metallic phase 
(fi > p c ) from the insulating phase (fj, < (i c ). 

symmetry topological critical 



class 


invariant C exponent v 


A 


-1 


1.18 ±0.05 


A 





1.21 ±0.03 


All 


-1 


0.99 ±0.02 


All 





1.02 ± 0.04 



TABLE I: Critical exponent v for the metal-insulator transi- 
tions indicated by circles in Figs. [1] and [3] The value C of the 
topological invariant at the insulating side of the transition is 
indicated, as well as the symmetry class: class A (quantum 
Hall effect) or class All (quantum spin-Hall effect). 



To characterize the metal-insulator transition we have 
performed a finite-time scaling analysis, along the lines 
of Refs. [H, H3- One- parameter-scaling in 4D requires 
that the mean-square-displacement obeys the scaling law 



A 2 (*) 

Vt 



= 2?(r 4 i), eocl/x-Mcl" 



(13) 



All microscopic parameters enter only through the local- 
ization length £, which has a power law divergence at the 
critical point \i c . We determine the critical exponent v 
by expanding I? in a power series near fj, c , including also 
finite-time corrections to scaling [28[. A typical scaling 
plot is shown in Fig. [2] 29] . 

Results are listed in Table U (first two rows). As ex- 
pected from general considerations there is no de- 
pendence of v on the value of the topological invariant C 
at the insulating side of the transition. Within numeri- 
cal accuracy, v w 1.2 for the transition from a metal to 
either a trivial or a nontrivial quantum Hall insulator. 

The calculations described so far are for a system with 
broken time-reversal symmetry, in class A of the quan- 
tum Hall effect. To investigate the role played by this 
symmetry we have repeated the calculations in the time- 
reversally-symmetric class All of the quantum spin-Hall 
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FIG. 3: Phase diagram of the 4D quantum spin-Hall effect 
(preserved time-reversal symmetry, class All). The colour 
scale is as in Fig. [T] 



effect, by replacing the t z matrix in Eq. © with the unit 
matrix ro. The potential V{x) then commutes with T 
and preserves time-reversal symmetry. The phase dia- 
gram is qualitatively the same, see Fig. [31 at least for not 
too strong disorder. The insulating phase for 1/K > 10 
is now connected to the C — phase at weaker disorder, 
so we definitely know that this is a trivial insulator. 

Tabled] (last two rows) shows that a significant effect of 
time-reversal symmetry appears in the critical exponent 
v w 1.0 for class All — well below what we found for 
class A. Within numerical accuracy, the critical exponent 
which we find for the 4D quantum spin-Hall effect agrees 
with the semiclassical formula UM Lis 13011 



d-2 



(14) 



for (i-dimensional Anderson transitions with time- 
reversal symmetry. 

In conclusion, we have demonstrated that the phase di- 
agram of the quantum Hall effect depends sensitively on 
the dimensionality. While in 2D only insulating phases 
exist, in 4D a metallic phase appears as well, for an in- 
termediate range of not too strong and not too weak dis- 
order. Topological phase transitions in the 4D quantum 
Hall effect occur via consecutive insulator-to-metal-to- 
insulator transitions, rather than via a single insulator- 
to- insulator transition as in 2D. This altogether different 
phenomenology is made possible by the fact that bro- 
ken time-reversal symmetry prevents the formation of a 
metal in 2D but not in higher dimensions. 

One direction for further theoretical research is to- 
wards a semiclassical scaling theory of 4D quantum Hall 
transitions. In 2D there is no hope for such a theory, 
but in 4D the situation is more promising in view of the 
close agreement between numerics and semiclassics that 
we have found for the critical exponent in the quantum 
spin-Hall effect. From the experimental point of view, the 
mapping of a scalar 3D Hamiltonian onto a ID dynami- 
cal system has been realized in cold atoms with a pulsed 
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optical lattice, driven by three incommensurate frequen- 
cies Our 4D-to-lD mapping adds the complexity 
of an additional driving frequency, but more importantly 
would require control over an internal degree of freedom 
of the atoms to produce the required 4x4 matrix struc- 
ture. 
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Appendix A: Finite-time scaling analysis of the 
metal-insulator transitions 

Here we give the details of the scaling analysis pre- 
sented in the main text, following the established proce- 
dure [28| , Close to the transition point /j, c , the anoma- 
lous diffusion coefficient V = A 2 /\/i obeys the single- 
parameter scaling law (|13j) . which we can write in the 
form 

\nV = F{t L ' iv u), u = n- fi c + 0((i- Ai c ) 2 . (Al) 

Finite-time corrections to scaling add a term 
t-yG^/^u) to In 2?, with y > the leading irrel- 
evant exponent. 

We expand the functions F and G in a power series in 
t 1 / 4,y u, and expand u in a power series in fi — fi Cl 

Nt N 3 

In V = lnV c + Y jC ^\t 1 ' iv u) k + g c^t^^u)*- 1 , 

fc=l k=l 

(A2) 

N 2 

u = f i- f i c + Y / c ( k\v-Vc) k - (A3) 

fe=2 

We fit the unknown parameters to the numerical data 
for A 2 (/i, t), averaged over 500 — 1000 initial conditions, 
varying /i around /i c for six values of t in the range 10 4 — 
6 • 10 7 . The number of fit parameters Ni,N2,Na was 
increased until we reached a value for chi-squared per 
degree of freedom (x 2 /ndf) close to unity. 

We have performed this scaling analysis for two sets of 
incommensurate frequencies, in order to make sure that 
our results for the critical exponents do not depend on 
the precise choice of driving frequencies. The two sets 
are 
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class C freq. set fi c Ni N2 N3 v x 2 / nc h" • Set fii: W2 = 2n/X, 0J3 = 2ir/\ 2 , u>4 — 2n/y/2, where 

A ~J f2 3 52 1 2 1 18 ± 05 97 ^ ^ s * ne rea ^ solution to the equation x 3 — x — 1 = [HI ■ 

A -1 fi 2 3J51 2 1 L21±(L06 L01 • Set fi 2 : = 2jln3, o; 3 = 27rln4, o; 4 = 27rlii5. 

Results tor both frequency sets are given m lableslllland 
IIIIl and summarized for frequency set Q\ in Tabic HI 



L2 

A Qi 4.30 1 1 1 1.21 ± 0.03 0.89 
A fi 2 4.30 1 2 1 1.20 ±0.03 1.05 



TABLE II: Parameters for the scaling analysis of the metal- 
insulator transition in the 4D quantum Hall effect (symmetry 
class A). 



class C freq. set fi c Ni N2 N3 



X 2 /ndf 



All -1 fii 3.65 1 1 1 0.99 ±0.02 1.16 

All -1 fi 2 3.62 111 0.98 ±0.04 0.78 

All fii 4.22 1 1 1 1.02 ± 0.04 0.88 

All fi 2 4.24 1 2 2 1.07 ±0.05 0.91 



TABLE III: Parameters for the scaling analysis of the metal- 
insulator transition in the 4D quantum spin-Hall effect (sym- 
metry class All). 



